Interface-mediated pairing in field effect devices 
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We consider the pairing induced in a strictly 2D electron gas (2DEG) by a proximate insulating 
film with polarizable localized excitations. Within a model of interacting 2D electrons and localized 
two-level systems, we calculate the critical temperature Tc as a function of applied voltage and for 
different materials properties. Assuming that a sufficient carrier density can be induced in a field- 
gated device, we argue that superconductivity may be observable in such systems. Tc is found to be 
a nonmonotonic function of both electric field and the excitation energy of the two-level systems. 
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Shortly after the publication of the BCS theory of su- 
perconductivity, W. A. Little proposed a possible pair- 
ing mechanism for electrons in long organic molecules in- 
volving localized electronic excitations in the molecules' 
side chains Jj. The suggested advantages of such an ar- 
rangement included the large characteristic energy of the 
couplings, which might lead to high temperature super- 
conductivity (HTSC), and the separation of the excita- 
tions themselves from the screening effects of the electron 
gas. Such a pair mechanism has never been realized, 
presumably due to large fluctuation effects in quasi- ID 
systems. Later, Ginzburg ^^id AUendcr, Bray and 
Bardeen 3] proposed a similar excitonic 2D mechanism 
in superconductor-semiconductor sandwiches. The enor- 
mous body of early work on this problem has been re- 
viewed in ^] . One problem with these schemes is clearly 
that the localized nature of the excitations implies that 
at most one or two atomic layers of the intercalating in- 
sulator can contribute to the pairing, meaning that the 
enhancement of pairing in the relatively thick nearby su- 
perconductor is negligible. 

Here we propose that a similar scheme might work for 
an insulating layer in proximity with a superconducting 
layer of near- atomic thickness. Such systems can in prin- 
ciple be prepared in several ways, but the most promising 
is perhaps the field effect method pioneered in the early 
90 's with the intent of increasing the carrier density and 
hence in HTSC cuprate devices [1, i, 0, i, II- This 
work has shown incontrovertibly that Tc and other super- 
conducting properties in metallic samples can be influ- 
enced by an applied gate voltage. While the experiments 
have most often been interpreted in terms of electronic 
structure variations near the interface, they are not un- 
derstood in detail. 

It is interesting to ask if new physics can result from 
field doping of insulators. In the most recent reports on 
field-doping of SrTiOg, Pallecchi et al. ^ and Ueno et 
al. pj| were able to achieve an areal carrier density of 
^ 0.01-0.05 per unit cell. It is difficult to achieve higher 
densities due to electrical breakdown in the insulating 
layer at high fields. But there seems to be no fundamen- 



tal objection to even higher charge densities, since com- 
plex oxide dielectrics and ferroelectric oxides can achieve 
polarizations in the range of ~ 0.5 Q. In fact, very re- 
cently the first observation of field-induced superconduc- 
tivity was reported for a device with a Ndi.^Bai.gCuaOy 
epitaxial film grown on SrTiOs substrates [12| . 

Therefore one may legitimately ask the question what 
critical temperature can be achieved in a system where 
the pairing comes entirely from the excitations in the 
proximate insulating layer, and how Tc is likely to vary 
with field in this case. Alternatively, one can assume 
the existence of a 2D superconducting layer with pre- 
existing pair interaction and bare critical temperature 

, and ask by how much the presence of the insulat- 
ing layer enhances Tc. While in a strictly 2D system the 
Tc's referred to cannot correspond to true long-range or- 
der, the creation of a field-tuned superconducting state 
with algebraic order at finite temperatures would be of 
considerable interest and applicability in small devices. 

We begin by considering an insulating amorphous film 
LI in proximity to a correlated insulator L2 (drain-source 
(DS) channel of a field effect device), similar to Ref . [l^. 
with a small density of localized charge carriers. Ap- 
plying an electric field as shown in Fig. [T] sweeps charge 
carriers to the interface with the film LI where they ac- 
cumulate [l^ . The formation of the 2D band induce d by 
an electric field has been studied by Poilblanc et al. 'hS\. 
We assume for the moment that in the absence of the film 
LI there are no pairing interactions between the electrons 
in the material L2. Qualitatively, we expect the follow- 
ing picture to apply: virtual excitations in the dielectric 
LI induce Cooper pairing in the adjacent layer L2. The 
critical temperature must increase initially with the field 
since carriers are being injected into the system. With 
increasing electric field the larger level splitting of the 
two level system leads to a suppression of the polariza- 
tion fluctuations and the pair potential decreases. Tc is 
therefore expected to reach a maximum at a characteris- 
tic field strength; it is our objective to estimate the scale 
of possible Tc's through this process, as well as the field 
strength required to attain it. 
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FIG. 1; Our model of a field efltect transistor is represented 
by two layers (LI and L2): 

LI: dielectric layer with local dipoles of dipole moment edsp, 

presented by two-level systems; 

1/2: metallic DS-channel with a 2D electron gas. 



We now propose a crude but concrete framework 
within which one can calculate these effects. Roughly 
speaking, the dielectric can have two effects on Tc. First, 
it can reduce the Coulomb pseudopotential of electrons in 
the field-doped layer due to the large dielectric constant 
of the amorphous insulator LI. We are more concerned 
here with the second effect, namely an additional contri- 
bution to the residual pairing interaction at the interface 
due to virtual polarization of the dielectric itself. 



I. MODEL 

The Hamiltonian we consider describes a single layer 
L2 of electrons cj hopping on a square lattice with lat- 
tice constant a, nearest neighbor hopping t, and a set of 
localized two-level systems in an adjacent layer LI with 
ground state sj and excited state pj separated in en- 
ergy by Asp, and their respective dipole moment is edgp 
(Fig. [T|) . Note these states simply designate ground and 
excited states of a localized two level system, and need 
not correspond to actual atomic s- and p- orbitals. To 
ensure that the two-level system is occupied by only one 
electron, the constraint slsi+plpi = 1 must be enforced. 
The electric field 



(1) 



is assumed to populate the metallic layer and simultane- 
ously polarize the two level systems by driving transitions 
between the s and p states, which are coupled to the free 
electrons in the layer L2 by a contact interaction Vsp. 

The Hamiltonian contains, according to our assump- 
tions, the following elementary processes: 
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Here Ht describes a band of noninteracting 2D electrons 
on a square lattice, H21 the energies of the localized two 
level system, H^xt the coupling of the electric field to 
these orbitals. Hint the Coloumb interaction between 
electrons in the metallic layer and the two level sys- 
tem, and the chemical potential. The direct electron- 
electron interaction term i/e-e is taken to be local and 
repulsive. Very similar models have been used recently 
to discuss dielectric properties of bulk cuprates [1^1 , the 
competition between charge density wave and supercon- 
ductivity in two dimensional electronic systems [iq , and 
superconductor-ferroelectric multilayers [l7j . The sys- 
tem is not exactly soluble, but will be treated under the 
assumption that the polarization through the charges in 
L2 is not large enough to drive the system into the ferro- 
electric state. We first diagonalize H21 + Hext, and then 
express the corresponding quasiparticle operators in a 
pseudospin representation where Sf measures the occu- 
pation of the 2-level system and Sf^ induces transitions 
between the two eigenlevels. The occupation constraint 
on the two-level system is preserved by the usual spin 
algebra. The (exact) final form of the Hamiltonian is 



H21 + Hext 
Hint 



-2^EI + {\AspYY.S'^ 

i 



(9) 
(10) 
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where 



K = 2K 



E. 
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El + {\Aspf 
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Ht + H2 



He 



Hint + Hn + He-e 



(2) 



The system of interacting spins- 1/2 and fermions may 
now be treated within linear spin wave theory, an ap- 
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proximation which is justified a posteriori by the ob- 
servation that the occupation of the higher-energy 2- 
level state is very small, due to the sufficiently large 
2-level splitting [isll . Introducing the usual Holstein- 
Primakoff bosons [l9| , we make the approximate replace- 



by rewriting the fourth term of Eq. (|T2 



ments 



b^jbj to find 



<ij>,a i.a i 

+ K cl^C^^^blbi + V^Y^ claC^,<T{b^ + b}) + H^-e 

(12) 



We have introduced the energy splitting of the two-level 
system 



E2 



El + (-A,p)2 



(13) 



and the chemical potential has been renormalized /i^ = 
/i + K/2. 

We now declare that the Hamiltonian H is as good a 
starting point as ©-([HI) provided <C 1. We there- 

fore proceed to apply a Feynman variational procedure 
to the exact i?, Eq. (fT^ . after a unitary transforma- 
tion, and attempt to extract the pairing mechanism. The 
dominant term for this mechanism is the term with Vx 
which, as in the phonon-induced superconductivity, pro- 
duces Cooper pairing in second order perturbation the- 
ory [if. However, for the field strengths which we will 
consider, the term with I4 is also not negligible and it 
will alter pairing in one significant respect: since it mod- 
ifies the 2-level splitting when a charge carrier occupies 
the site coupled to the considered 2-level system, it will 
either increase the pairing interaction (for reduced split- 
ting, i. e., negative Vz) or decrease pairing (for enhanced 
splitting, i. e., positive Vz) — an observation which may 
be derived straightforwardly from a mean field decou- 
pling. For the physical system that we consider, E^p and 
Vsp always have the same sign, and consequently Vz is 
positive [2Q|. All of the above considerations can also 
be derived in terms of states where charges in LI are lo- 
calized in states of definite position with respect to the 
interface, as depicted schematically in Fig. 1. In this 
representation, the Vz interaction may be shown to cor- 
respond to the repulsion of the electron in layer L2 and 
the field- induced dipole in LI. 

In order to incorporate this effect we will include below 
the thermal average VzY.^,^{c\^^Ci^a)b\bi = VznY,ib\bi 
explicitly in the energy of the two-level system before 
handling the bosonic degrees of freedom. This is achieved 



aC^,<jblb, 



= -NVzUrib + VzU Y ^Ibr + VzTlb ^ 4<t^ 



Vz Y ~ n/2 



rib 



(14) 



where 71^ — {blbi) is the number of bosons (the rela- 
tive number of inverted 2-level systems per site), n = 

(l/N) J2i ai^icr^i,'^) ^hc charge carrier density in layer 
L2, and N is the number of sites in L2. The second term 
on the right hand side accounts for a density-dependent 
renormalization of the two-level splitting: 



^21 



E, 



nVz 



(15) 



In a first step we apply a Lang-Firsov (LF) transfor- 
mation of the Hamiltonian H in order to identify the 
pairing interaction, and in a second step we control the 
decoupling of the interaction terms through Feynman's 
variational principle which also fixes the parameters of 
the LF transformation. 

The LF transformation of the Hamiltonian H = U^HU 
is achieved by the following unitary operator U: 



with 



f/ = exph5i(0)] exph52(7)] 



Si{9) = -^^E(^l-^^ 



i.rr 



(16) 

(17) 
(18) 



The parameters 9 and 7 will be fixed through the vari- 
ation of the free energy. The standard form of the 
LF transformation for the Holstein model takes = 
and 7 = 1. With a "zero-phonon" approximation the 
(phononic) Holstein model accounts for exact results in 
the antiadiabatic limit E2i/t 3> 1. The variation of 
the parameters 9 and 7 has been devised in order to 
reproduce the adiabatic limit of the Holstein model as 
well [21I, m, m, ll^l ; the static case is realized with 7 — > 
and with a finite 9 which signifies a displacement field. 
For the Holstein model, a second canonical transforma- 
tion of Bogohubov type Ub = exp —aj^ii^l^l + ^j^O 
is sometimes applied to allow for the anharmonicity of 
the lattice fluctuations. For our model, we have to sup- 
press states with two bosons (corresponding to our initial 
constraint for the 2-level system), and subsequent vari- 
ation of the free energy shows that in fact a = 0. Also 
the "static displacement" should be suppressed since we 
do not intend to consider a finite static polarization. We 
verified from the minimization that indeed ~ for the 
considered range of small to intermediate T4p/4i. For 
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this reason, and to simplify the notation, we fix 9 and a 
to zero. 

With the LF transformation U — exp [—£'2(7)] we ob- 
tain 



H^Eq + Hi^f+ Hv^ + Hv^ + Hv^v^ (19) 
where — —^NE2i — NVzUni,, and 

-lib}-bj) 



e 21 e 



<iJ>,CT 
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i,(T i 



with 



(1 - 7) K 



H 



■i.a 



(20) 

(21) 
(22) 



z \ ^ 



^2/ ■ 



(23) 



where the effective chemical potential is given by 



/^LF=/i + VJ2 + 7(2-7)^ 

^2/ 



14^6 (24) 



The induced electronic interaction Vcs is now 



y2 

Fcff = 2-^7 



2/ 



(2-7)-7(3-«)^ 



U , 



(25) 



where we note that the Hubbard term i?e-e has been 
unaffected by the above transformations of the electron- 
two level system coupling terms, due to the fact that it is 
simply a product of local densities. The local repulsion 
therefore merely diminishes the effective attraction by U . 

If 7 = 1, the first term within brackets in Eq. ([25]) coin- 
cides with the well-known result from the exact mapping 
of the Holstein model to the attractive Hubbard model in 
the high-frequency limit {E21 — > 00). The second term in 
the brackets oc Vz is always repulsive and suppresses Tc 
for increasing field. The electric field dependence enters 
through Vx, Vz, E21 and the implicit dependence of 7 on 
Esp . 

We have now achieved the desired result of expressing 
the Hamiltonian in terms of an effective BCS-like interac- 
tion between electrons, at the cost of introducing bosonic 
phase factors into the hopping matrix elements. 
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FIG. 2: Left panel: Effective pair interaction Vcs/At vs. nor- 
malized electric field energy Esp/4:t for various excitation en- 
ergies Asp of the local two- level systems in layer LI. The 
dielectric constant e — 100 controls the increase of charge den- 
sity in the metallic layer L2 with electric field (cf. Eqs. (|38p - 
(|39p ). whereby we chose a bandwidth 4t = 400 meV, a dipole 
length dsp = 2 A (for charge transfer excitations), and a lat- 
tice constant a = 4 A. The coupling of the dipoles to the 
conduction electrons is fixed at Vsp/4:t = 1.89 which corre- 
sponds to a spatial distance r/a — 1.5 (cf. Eq. (|40p ). Right 
panel: variational parameter 7 vs. normalized electric field 
Esp /it, where 7 is found from the minimization of the free 
energy (rhs of Eq. ((26])). 



We begin by setting [7 = and investigate the mag- 
nitude of the attractive interaction obtainable by polar- 
izing the two level systems. In Fig. [2] we illustrate how 
Vcs depends on applied electric field. The main point is 
physically obvious: if the electric field is sufficiently large, 
it polarizes the electric dipoles in LI and suppresses the 
Little-type mechanism. We also note that if the bare ex- 
citation energy A^p of the dipoles in the insulating layer 
is small, the intrinsic low-field pairing strength can be 
quite large, up to several electron volts. However Asp 
should not be set to considerably smaller values than in 
Fig. [2] in order to guarantee the condition that the occu- 
pation of the excited state is negligible. 

In the second step we introduce an exactly solvable 
test Hamiltonian i/tcst and determine the fields in TJtcst 
through a Bogoliubov inequality for the free energy F of 
the model system [25j |: 



F < Fu 



{H — i?tcst)t 



(26) 



where ( )tost signifies the thermodynamic average with 
the test Hamiltonian. We already anticipate the result 
of this variational scheme |29| and write 



E, 



O,tost 
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(27) 



5 



0.6 



0.4 



0.2 



.... A /4t = 


3.75 


sp 




A / 4t = 


2.5 


sp 




__ A /4t = 


1.25 


sp 






le-06 



le-08 



le-10 = 



le-12 



- le-14 



0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

E I At 



FIG. 3: Left panel: Effective hopping tcs/t vs. normalized 
electric field energy Eap/it at the transition temperature Tc; 
all paramaters are identical to those of Fig. [2] Right panel: 
bosonic occupation number nt, (the relative number of in- 
verted two-level systems per site) vs. normalized electric field 
energy Esp/At at the transistion temperature Tc. 



with 



N 2 

-Ectest — Eq — —n VcS 



t exp 



(28) 
(29) 



We have assumed s-wave pairing for simplicity. The as- 
sociated gap equation then reads 

— — - tanh — ^ 
2Ek 2T 



A = Foff ^ TTT^ tanh ■ 



(30) 



and the Bogoliubov quasiparticle dispersion 



(31) 



—2 toff (cos kx + COS ky) 

To understand qualitatively the physics of the corre- 
lations induced by interaction of the metallic layer with 
the 2-level systems in LI, we display in the left panel of 
Fig. [3] the magnitude of the renormalized hopping tcff • 
The band narrowing can be rather significant for small 
excitation energies A^p but for the parameters of greatest 
interest will turn out to be only a factor of 1-5. 

The fields still depend on the parameter 7; correspond- 
ingly the rhs of relation Eq. (|26p has to be minimized with 
respect to 7. For this purpose we need the explicit form 
of Ftost: 

1 



Fto.st — El 



O.tcst 



^ iV In (1 - e 



with 



BCS 



f3 



^In, 
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BCS 



(32) 



5l(^k-a + A(c_fe,4Cfe,T)) (33) 



Since we have chosen i?o,tcst so as to guarantee the re- 
lation (iJ — i?tcst)tcst = 0, the variation reduces to finding 
the optimal value of 7 from the minimization of Ftest- 



n{n + 2) 



n{n + 2) + n{n + 1){2 -n){V J E*i) + S{j) 



where 



5(7) 



^2/ 



coth(^) 



and 



St 



£k 



(34) 



(35) 



(36) 



For the range of parameters considered in our evaluation 
the function 5(7) can be neglected with respect to the 
other terms in the denominator of Eq. . This allows 
us to calculate 7 algebraically from Eq. ((34)) . The valid- 
ity of this approximation has been proved by iterating 
the implicit Eq. Fig. [2] (right panel) displays the 

decrease of 7 with increasing Esp for various values of 
A^p at fixed Vsp/U = 1.89. 

Finally the assumption below Eq. pT|) that the bosonic 
occupation number ni, — {b\bi) is very small is in fact 
justified a posteriori for temperatures up to the transition 
temperature Tc. In the considered parameter range of the 
two- level splitting Agp and of the applied gate field E'sp, 
rifc is always smaller than 10~^. Typical field dependences 
of Tif, are displayed in the right panel of Fig. [3l 
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FIG. 4: The derivative of the chemical potential with respect 
to particle density, dn/dn, versus particle density n at the re- 
spective Tc, calculated in Sec. II. All paramaters are identical 
to those of Fig. [2] For a positive derivative, the normal and 
superconducting states are globally stable in the vicinity of 
the transition. 

The effective model does not implement the long-range 
part of the Coulomb interaction. A model with purely lo- 
cal interactions, however, has a tendency towards phase 
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separation. In order to investigate if the transition from 
the normal to the superconducting state is a transition 
between global minima of the free energy, we evaluated 
the derivative of the chemical potential with respect to 
carrier density. In this work, we focus on the transi- 
tion into the superconducting state and not on the evolu- 
tion of the superconducting state at lower temperatures. 
Then it is sufhcient to discuss the stability in the vicinity 
of Tc. Fig. |4] illustrates our conclusion that the derivative 
of the chemical potential with respect to particle density 
is always positive for A^p/At larger than approximately 
1.4. Correspondingly, the transition into the supercon- 
ducting state is not preempted by a competing transition 
into a phase separated state. At somewhat lower values of 
Asp/4t, the phase separation will probably be suppressed 
by the non-local part of the Coulomb interaction. As 
discussed in the following section, the relevant values of 
Asp/4i for a sufficiently strong interface-mediated pair- 
ing appear in the regime where djji/dn is positive. 



II. RESULTS 

In order to present the numerical results in the phys- 
ically relevant range of control parameters, we have to 
relate microscopic quantities to laboratory parameters, 
such as the dielectric constant e of the insulating layer, 
the dipole length dsp of the two-level systems, and the 
distance r from the two-level systems to the sites of the 
conducting layer. Furthermore we have to establish a re- 
lation between the electric field energy E^p (see Eq. ([T])) 
and the number of charge carriers in L2. In our calcula- 
tion we assume that the charge carrier density is a linear 
function of the electric field with a field independent ca- 
pacitance C of the dielectric LI: Q — CV, where V is 
the voltage drop across the dielectric and Q is the total, 
accumulated charge at the interface in L2. Then, the 
charge per square unit cell is 



n^C 

e A 



(37) 



where a is the lattice constant, A is the area of L2 and d 
its thickness (cf. Fig. [T]). With C = CocA/d and Eq. ^ 
we establish 



n = c 



IT 



with 



6 dsv 



(38) 



(39) 



The interaction energy between the dipoles next to the 
interface and the electric field of a charge carrier on the 
nearest site is 



1 e^d,, 
47reo 



We note that dynamical screening of the interaction 
Vsp should be included in order to get more precise es- 
timates. Although this will reduce somewhat the high 
values of Tc in our evaluation, it will not alter the qual- 
itative behavior of the Tc-dependcnce on the gate field 
(see below). 

Finally, the excitation energy of the dipoles has to be 
identified. A generic dielectric is not composed of a single 
species of well-defined 2-level systems. Dipole excitations 
at various energies are always present. We now address 
excitations in three different energy ranges: for 10 eV 
(low atomic excitations with dipole length dsp ~ 1 A), 
in the 1 eV range (charge transfer excitations with d^p ~ 
2 A), and in the meV range (ionic displacement in atomic 



clusters with d 



sp 



0.1 A). 



(40) 



A. Ionic displacement 

In dielectrics a displacement of ions, well localized at 
atomic positions, usually accounts for the high dielec- 
tric constant. Such displacements in atomic clusters 
typically correspond to a dipole length of the order of 
dsp ~ 0.1 A and a small excitation energy. In this case, 
the excitation energy is related to e and dsp through 
Asp = 87re^fisp/(a^ (e - 1)). This relation follows di- 
rectly from the polarization density (P) of the dielectric 
(with a cubic unit cell of volume a^) in linear response: 
(P) = edsp{slp^+pls,)/a^ = {2 e'^d%/a^ Asp) S where 
(P) = (e — l)/47rf holds. Ionic displacements may have 
energies of the order of 10 meV (for e of the order of 20). 
However, the induced pair interaction is far too small 
(see Fig. [5]) to find sizable transition temperatures. It is 
the repulsive term in VeS / it of Eq. (|25p which very effec- 
tively impedes a transition to the superconducting state 
for such small values of Agr,. 



B. Charge transfer excitations 

In order to realize a sizable Tc of the order of 100 K, 
excitations at intermediate energy and sufficiently large 
dipole length have to be available for polarization. 
Charge transfer excitations in the dielectrics are found 
at these energies — e.g. for the transition metal (TM) 
oxides at energies of the order of 1 eV, the charge trans- 
fer gap. Although these excitations of, for example, a 
TM-oxygen plaquette or octahedron become delocalized 
through hybridization, we assume that most states near 
the interface have been localized by disorder and local 
strain 27]. Furthermore localized bound states below 
the charge transfer gap are conceivable. It is not the 
focus of our present investigation to identify these local 
excitations for a specific material. We rather point out 
that such processes are possible and then evaluate Tc as 
a function of their respective energy. 

We intend to focus on this latter type of localized 
charge transfer excitations. For the evaluation we now 
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FIG. 5: Ionic displacement in atomic clusters: induced pair- 
ing interaction Ves/it versus electric field energy Esp/4t for 
e = 10, dsp = 0.1 A and r/a — 1 (which corresponds to 
Vsp/4:t — 0.2125). Left and lower right panels show different 
scales. Upper right panel: Pairing interaction versus dielec- 
tric constant e for charge density n = 0.5. 



fix the following parameters: the bandwidth, the dipole 
length, and the lattice constant: 
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FIG. 6: Charge transfer excitations: transition temperature 
Tc/At versus excitation energy Asp/4/:. The three curves 
present Tc for small and intermediate band filling. The di- 
electric constant is e = 100, which implies c = 1.87 for the 
considered set of parameters (see Eq. (|4T])). The coupling of 
the dipoles to the conduction electrons is fixed to Vap/'U = 
1.89 which corresponds to a spatial distance r/a = 1.5 (cf. 
Eq. (go])). 



At = 400 mcV 
dsp — 2 A 
a = 4A 



(41) 



Here, we chose a bandwidth and lattice constant typi- 
cal of the high- Tc cuprates and many other oxides, and 
we take a dipole length which corresponds to the inter- 
atomic oxygen-TM distance. For excitation energies A^p 
of more than 0.2 eV, the primary polarization processes 
are electronic in nature. 

In Fig. [6] the transition temperature to the supercon- 
ducting state is displayed as a function of the excita- 
tion energy. In the range of small, increasing values of 
Asp we observe a strong enhancement of Tc whereas the 
transition temperature decreases with excitation energies 
Asp/4t above « 2.5. The latter observation is expected 
since the pairing interaction VcS is inversely proportional 
to the excitation energy for large Asp. For small Asp, the 
repulsive term (second term in Vcs-, cf. Eq. ([^5]) ) domi- 
nates and suppresses the transition to the superconduct- 
ing state for a finite value of Asp. Note that the finite 
longitudinal pseudospin-charge Vz is decisive for the de- 
cay of Tc at small excitation energies, outside the regime 
where the Holstein model is appropriate. 

For small increasing electric field, the transition tem- 
perature is raised due to the accumulation of charge in 
the metallic layer (cf. Fig. [7]). The electric field strength is 
directly related to the band filling or induced areal charge 
density. Strong electric fields with sizable band filling 
lower the transition temperature as the repulsive term in 
14ff is enhanced and, moreover, the effective level split- 
ting is enlarged. In fact, as seen from the Asp/4t — 1.25 
curve in Fig. [7] there are two scales for the suppression of 
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FIG. 7: Charge transfer excitations: transition tempera- 
ture Tc/At versus electric field energy Esp/At and band fill- 
ing n, lower axis and upper axis, respectively. The three 
curves present different dielectrics which are characterized 
by different excitation energies only, other parameters are 
fixed for comparison. The dielectric constant e is 100 and 
Vsp/U = 1.89 (r/a = 1.5). 



Tc at higher fields. The lower scale is set by the repulsive 
term in Vea and is also responsible for the observed decay 
of Tc in the two further curves (with Asp/4t = 2.5 and 
3.75). In this regime, Tc is being suppressed primarily by 
the increasing repulsion Vz between the polarized dipoles 
and the 2D electrons, which scales with the applied field. 
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FIG. 8: Charge transfer excitations: transition temperature 
Tc/At versus excitation energy Asp/4i. The three curves 
present diflterent dielectrics (e = 10,50, 100). The band fill- 



ing is fixed at n = 
corresponds to r/a 



0.1, and the interaction 
= 1.5. 



Vay/U = 1.89 



The larger scale, which is responsible for the slow decay 
at even higher fields, is set by Agp/it, and corresponds 
to the eventual saturation of the dipole moment of the 
2-level systems. 

The nonmonotonic dependence on filling or electric 
field is reflected in the varying height of the three dif- 
ferent curves in Fig. [HI However more striking is the 
small variation of the position of the maxima in Fig. [6] 
with fining (from 0.01 to 0.5). A value of A^p ~ 1 eV 
seems to be optimal for a bandwidth of At = 400 meV. 
This optimal excitation energy in units of the band width 
is approximately 



A°PV4i -2.5 



(42) 



A°P* is weakly dependent on the parameters Vsp, dsp 
and e which mostly influence the maximum value of T^. 
For fixed charge density, dielectrics with larger e display a 
higher transition temperature (for given A^p), see Fig.|Sl 
as the electric field necessary to create the charge density 
is smaller, accounting for a larger Veff- 

An increasing dipole charge-carrier interaction Vsp is 
not only responsible for an enhancement of Tc but, for 
sufficiently large Vsp, also for a "retarded" initial increase 
of Tc with electric field (see Fig. [S]). Again, this obser- 
vation may be traced back to the field and interaction 
dependence of the repulsive term in Ves. 

For direct comparison with experimental devices, we 
note that the gate voltage is related to the electric field 
energy plotted in the figures through Esp/At — edspV/d 
where d is the thickness of the dielectric gate. As a con- 
crete example for a dielectric gate layer, we take SrTiOa. 
Breakdown fields of 4x10'' V/m are reported for this ma- 
terial with a low-T dielectric function of order 100 [1^. 
The maximum of the Tc versus electric field energy curves 
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FIG. 9: Charge transfer excitations {dsp = 2 A, Asp/4i = 2.5, 
e = 100): transition temperature Tc/4t versus electric field 
energy Esp/4t, parameterized by r/a. The corresponding Vsp 
is in the intermediate coupling range: Vsp/4:t — 4.25 (corre- 
sponds to r/a = 1), Vsp/4t = 2.72 (r/a = 1.25), Vsp/At = 1.89 
[r/a = 1.5), Vsp/4:t = 1.39 (r/a = 1.75). 



occurs for n = 0.2 at about 2x10* V/m. Thus the re- 
quired fields for the maximum Tc are only about 5 times 
higher than those already realized in this system. Given 
that breakdown occurs due to "pinhole" -type defects in 
the films, it seems to us that the manufacture of samples 
with the required breakdown fields is challenging but far 
from impossible. 

We conclude this subsection with a brief discussion 
of the consequences that a nonzero local interaction U 
within the charge-carrier layer L2 has on the reduction 
of Tc in a weak-coupling evaluation. As shown in Eq. 
(PS)) . such a repulsive Hubbard- type interaction docs not 
modify the field dependence of the effective interaction 
VcS_, but just adds a constant —U. For weak interac- 
tion U/At < 1, the effective attraction Vcs is reduced 
but not fully suppressed for small to intermediate fields 
(cf. Fig. [2]) and corresponding filling. This observation 
is refiected in the field and filling dependence of Tc (see 
Fig. [TU|) : Tc is reduced for small to intermediate fields 
and suppressed for strong field-induced doping. As a 
consequence of the field-dependent reduction of Tc, the 
maxima are shifted towards lower values of filling, from 
about n = 0.2 at [/ = to n = 0.1 at U/4t = 1. The 
required fields for the maximum Tc are reduced similarly. 



C. Atomic excitations 

At high energies, we encounter atomic excitations, such 
as the 2p to 3s transition of states in the metal ox- 
ides. These states are treated in exactly the same way in 
the current theory, but are characterized by small dipole 
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FIG. 10: Charge transfer excitations: transition temperature 
Tc/'it versus electric field energy Esp/At and band filling n, 
lower axis and upper axis, respectively. A local electronic 
interaction U (within layer L2) is included in a weak coupling 
evaluation. The dielectric constant e is 100 and Vsp/4f = 1.89 
(r/a = 1.5). 



FIG. 11: Atomic excitations (dsp = 1 A, Asp = 10 eV, e = 
100): transition temperature Tc/'it versus electric field energy 
Esp/4:t, parameterized by r/a: Vsp/'U — 2.125 (corresponds 
to r/a = 1), Vsp/4t = 1.36 (r/a = 1.25), Vsp/4t = 0.94 
(r/a= 1.5). 



lengths dsp ~ lA and large 2-level splittings Asp ^10 
eV. Polarizing such dipoles is difficult and T^s are cor- 
respondingly small. Comparison of Fig. [11] for parame- 
ters consistent with atomic polarizations with Fig. [5] for 
charge transfer excitations makes the distinction of the 
two scenarios evident. In the atomic case, Tc is so small 
that an observable effect may be found only for the small- 
est value of r/a. Moreover, the excitation energy is so 
high that a decrease of Tc is not seen, even for the largest 
electric fields. Of course, when band filling is above half 
filling, the transition temperature will decrease with in- 
creasing field. However these field strengths are beyond 
electrical breakdown. 



III. CONCLUSIONS 

We have proposed that the field-induced 2DEG in a 
field effect device may become superconducting entirely 
due to pair interactions with a proximate insulating layer 
with high polarizability. In the general case, properties 
of a DS-channel embedded in an oxide field effect tran- 
sistor will then be controlled not only by the electronic 
and structural properties of the DS-channel, as is usually 
assumed, but can furthermore be strongly influenced by 
the gate insulator or by other adjacent dielectric layers. 
The choice of the gate dielectric layer may then influ- 
ence the device behavior of a field effect transistor far 
beyond controlling the maximum gate polarization and 
gate current. 

For a concrete model of electrons in a 2D layer in- 
teracting with 2-level systems at the interface with the 



insulator, we calculated the superconducting critical tem- 
perature, which displays as a function of applied field a 
steep initial rise and subsequent decay. The rise is caused 
by the increasing density of charge carriers swept to the 
interface by the electric field and the decay is due to the 
interaction of the field-induced dipoles with the charge 
carriers. The optimal values of the field or gate voltage, 
as well as sample dimensions and dielectric properties, 
were discussed in some detail within the framework of 
this model. 

The goal of these model calculations was to establish 
the plausibility of superconductivity in field-effect devices 
enhanced or induced by the presence of a polarizable in- 
terface, and to investigate the magnitude of the criti- 
cal temperature and likely dependence on external field 
and materials properties. We recognize, however, that 
several potentially important aspects of the physics are 
not present in the model used. These include long-range 
Coulomb interactions, dynamics of the dielectric screen- 
ing, and correlations in the 2DEG. We anticipate that the 
effects of long-range Coulomb interactions are less impor- 
tant in the present case than in 3D metallic superconduc- 
tors due to the larger dielectric constants in the insulating 
material, but this must clearly be justified by a legitimate 
calculation of the true Coulomb pseudopotential entering 
the MacMillan formula. Screening of the bare electron- 
dipole interaction V^p must also be included; we expect, 
however, that the form of the dependence of Tc on the 
electric field will not change significantly (see Fig. [5]). 
The inclusion of a repulsive local interaction U does not 
qualitatively change this form either. However it reduces 
Tc (by a factor of 2 for U/U = 0.8) and it shifts the opti- 
mal value of the field-induced doping to lower levels (cf. 
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Fig.nni). 

If the effective on-site Coulomb repulsion in the drain- 
source layer is much larger than we have considered, s- 
wave superconductivity will be completely supressed. It 
is still tempting, however, to regard the starting Hamil- 
tonian ([3])-([7|) as the hopping of correlated electrons (as 
in the t — J model) where Coulomb interactions have al- 
ready been accounted for. In this case, however, one must 
implicitly assume that doubly occupied sites have been 
projected out, so that non-retarded s-wave superconduc- 
tivity is impossible. Modelling superconductivity within 
this framework in higher- angular momentum pair chan- 
nels will require including interactions among the two- 
level systems at the interface, so as to produce a nonlocal 
pair potential. While we have not calculated these ef- 
fects explicitly, as they are technically significantly more 
difficult, our expectation is that the creation and even- 
tual suppression of the paired state by the electric field 
will be qualitatively similar to what we have calculated. 
This expectation is supported by the observation that an 
increase in with doping has to result trivially from 
an enhanced carrier density. Moreover, the buildup of a 
field-dependent repulsive interaction Vz has to take place 
due to the eventual saturation of the dipole moments 
with electric field, even when the two-level systems inter- 
act. Investigations along these lines are in progress. 

In addition, our considerations can in principle be ap- 
plied to other, more weakly correlated systems than the 
copper oxides. Our model may be used to treat situations 
in which the DS-channel contains an attractive s-wave in- 



teraction besides the interface-mediated pairing. In this 
case we expect a qualitatively similar behavior with a 
nonmonotonic Tc- 

According to our calculations, devices with thickness 
of the insulating layer of order 1000 A subjected to 
voltages of order ^ 20V could display superconductiv- 
ity, if an optimal material can be found with a suffi- 
ciently large low frequency dielectric constant of order 
30-100 and strong quasilocalized electronic modes of en- 
ergy Asp/4t ^ 2.5 [11. With a 20 V gate field and a 
typical 18 nF capacitance of a 0.03 cm^ gate dielectric 
(1000 A), this would correspond to a surface charge den- 
sity of roughly 12 fiC /cm^. It is now routine to fabricate 
epitaxial high e oxide dielectrics, such as SrTiOa with po- 
larizations in the range of 10-40 /iC/cm^ ^8]. Interface- 
mediated 2D superconductivity therefore seems to us to 
be plausibly within reach if good interfaces can be man- 
ufactured. 
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